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Abstract. The topological analysis from Bjorkman (1995) for the standard model that describes the winds from hot stars 
by Castor, Abbott & Klein (1975) has been extended to include the effect of stellar rotation and changes in the ionization of 
the wind. The differential equation for the momentum of the wind is non-linear and transcendental for the velocity gradient. 
Due to this non-linearity the number of solutions that this equation possess is not known. After a change of variables and the 
introduction of a new physically meaningless independent variable, we manage to replace the non-linear momentum differential 
equation by a system of differential equations where all the derivatives are explicitely given. We then use this system of equations 
to study the topology of the rotating-CAK model. For the particular case when the wind is frozen in ionization (S = 0) only one 
physical solution is found, the standard CAK solution, with a X-type singular point. For the more general case (S t 0), besides 
the standard CAK singular point, we find a second singular point which is focal-type (or attractor). We find also, that the wind 
does not adopt the maximal mass-loss rate but almost the minimal. 
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1. Introduction 

Since the launch of the first satellite with a telescope on board, 
it has been established the widespread presence of stellar winds 
from hot stars. These winds are driven by the transfer of mo- 
mentum of the radiation field to the gas by scattering of radi- 
ation in spectral lines (Lucy and Solomon, 1970). The theory 
of radiation driven stellar winds is the standard tool to describe 
the observed properties of the winds from these stars. Castor, 
Abbott and Klein (1975, hereafter "CAK") obtained an analyt- 
ical hydrodynamic model for these winds, based in the Sobolev 
approximation. The CAK model has been improved by Friend 
and Abbott (1986,"FA") and Pauldrach et al. (1986,"PPK"), 
giving a general agreement with the observations. For a ex- 
tended review see Kudritzki and Puis (2000, "KP") and refer- 
ences therein. 

This agreement with the observations led to the devel- 
opment of a new method to determine galactic distances 
using Supergiants as targets, namely the Wind Momentum 
Luminosity relationship ("WML", Kudritzki et al. 1999, KP 
and references therein). 

More detailed studies from Puis et al. (1996) and Lamers 
& Leitherer (1993) came to the conclusion that the line-driven 
wind theory shows a systematic discrepancy with the observa- 
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tions. Lamers & Leitherer (1993) suggest that this discrepancy 
may arise due to an inadequate treatment of multiple scattering. 
Abbott & Lucy (1985), Puis (1987) and Gayley et al. (1995) 
have shown that multiple scattering can provide an enhance- 
ment of the wind momentum over that from single scatter- 
ing only by a factor of two - three for O stars (Abbott & 
Lucy (1985) found a factor of 3.3 for the wind of ^ Pup). 

Vink et al. (2000) calculate, including the multiple scatter- 
ing effects, mass-loss rates for a grid of wind models that cov- 
ers a wide range of stellar parameters. They found a much bet- 
ter agreement between theory and observation, concluding that 
the inclusion of multiple scattering increases the confidence of 
the WML relationship to derive extragalactic distances. In all 
the calculations involved in the WML relationship, the solution 
of the improved (or modified) CAK wind (hereafter m-CAK) 
is not used. Instead an ad-hoc yS-field velocity profile is utilized 
(see KP, Vink et al. 2000). 

The unsatisfactory results of the velocity field obtained 
from the m-CAK model when applied in the WML relation- 
ship could come from the complex structure of this non-linear 
transcendental equation for the velocity gradient and its solu- 
tion schema. Due to this non-linearity in the momentum dif- 
ferential equation, there exist many solution branches in the 
integration domain. A physical solution that describes the ob- 
served winds must start at the stellar photosphere, satisfying 
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certain boundary condition and reach infinity. There is no so- 
lution branch that covers the whole integration domain, thus a 
solution must pass through a singular point in order to match a 
second solution branch. Therefore, the solution in this second 
solution branch reaches infinity. To find the location of singular 
points is one of the most difficult aspects of topological analy- 
sis of non-linear differential equations. 

Bjorkman (1995) performed a topological analysis of the 
CAK differential equation. He showed that the original solu- 
tion from CAK, which passes through a X-type critical point 
and has a monotonically increasing velocity field, is the only 
physical solution that satisfies the condition of zero pressure 
at infinite radius. In this study Bjorkman did not include the 
influence of the star's rotational speed. 

Although it is known that the line-force parameters (see 
below) are not constant through in the wind (Abbott 1982), the 
standard m-CAK model still uses these parameters as constant. 
For the particular case of extreme low metalicities, Kudritzki 
(2002) introduced a new treatment of the line-force with depth 
dependent radiative force multipliers. As a test, he applied this 
new treatment for the most massive and most luminous O stars 
in the Galaxy and in the Magellanic Clouds (due to the lower 
metalicity) finding an acceptable agreement between theory 
and observations. Then, it was used to understand the influ- 
ence of stellar winds on the evolution of very massive stars in 
the early universe and on the interstellar medium in the early 
stages of galaxy formation. 

On the other hand, it is known from observations, that 
all early type stars have moderate to large rotational speeds 
(Hutchings et al. 1979, Abt et al. 2002) and for Oe and Be 
stars, their rotational speed is a large fraction of their break-up 
speed (Slettebak 1976, Chauville et al. 2001). The incorpora- 
tion of rotation in the CAK and m-CAK models has been stud- 
ied by Castor (1979), Marlborough & Zamir (1984), FA and 
PPK, concluding that the effect of the centrifugal force results 
in a downstream-shift of the position of the singular point, a 
slightly lower terminal velocity and a slightly larger mass loss 
rate as a result of an increasing in the star's rotational speed. 
Maeder (2001) studied the influence of the stellar rotation in 
the WML relationship, finding just a very small effect on it. 

A revision of the influence of the stellar rotation in radia- 
tion driven winds has been done by Cure (2004) finding that 
there exists a second singular point in these winds. He studied 
the case when the stellar rotational speed is high and found nu- 
merical solutions, that pass through this second singular point, 
and which are denser and slower than the standard m-CAK so- 
lution. 

In view of these results, it is crucial to understand the solu- 
tion topology of the standard model, forall when one wants to 
incorporate other physical processes into the theory. 

The purpose of this article is to study the topology of the 
rotating-CAK model. In section|2lwe give a brief exposition of 
the radiation driven winds theory and the non-Unear differen- 
tial equation for the momentum, including rotation, is shown. 
In section |3l after a coordinate transformation, we develop a 
general method that allows to replace the non-linear momen- 
tum equation in a simple and straightforward manner by a sys- 
tem of ordinary differential equations, where all the derivatives 



are explicitely given. In section ^ a general condition for the 
eigenvalue of the problem is developed. This condition allows 
to classify the topology of the singular point (Saddle or Focal) 
and constrains the location of it in the integration domain. 
Sectionlsjis devoted to the application of the criteria developed 
in sectionl^for the rotating-CAK model. In section|6lwe show 
numerical results of this topological analysis, first for a wind 
frozen in ionization (setting the 6 parameter of the line-force to 
zero) and compare our results with Bjorkman's (1995) for the 
non-rotating CAK model and the rotating-CAK model from 
Marlborough & Zamir (1984). Furthermore, section |6l shows 
the results of the influence of changes in the wind's ionization 
structure {6 + 0) on the topology and discuss the rotating-CAK 
wind model. Conclusion are summarized in section^ 

2. The non-linear differential equation 

The standard stationary model for radiation driven stellar winds 
treats an one-component isothermal radial flow, ignoring the 
influence of heat conduction, viscosity and magnetic fields (see 
e.g., Kudritzki et al. 1989, "KPPA"). 

For a star with mass M, radius 7?,, effective temperature T 
and luminosity L, the momentum equation with the inclusion 
of the centrifugal force due to star's rotation, reads: 

dv I dp GM(l-r) v2(r) 

V— = — + + g""'(p, v',nE) (1) 

dr p dr r 

where v is the fluid velocity, v' = dvjdr is the velocity gradient, 
p is the mass density, M is the star's mass loss rate, p is the 
fluid pressure, = VroiRt/r, where Vroi is the star's rotational 
speed at the equator, F represents the ratio of the radiative 
acceleration due to the continuous flux mean opacity, os, 
relative to the gravitational acceleration, i.e., F - o-gL/AncGM 
and the last term g''"*" represents the acceleration due to the 
lines. 

The standard parameterization of the line-force (Abbott, 
1982) reads: 

g''"' = ^ foir, V, v') (r^v v'f (nE/W(r))' (2) 

W(r) is the dilution factor and foir, v, v') is the finite disk cor- 
rection factor The line force parameters are: a, 6 and k (the 
last has been incorporated in the constant C), typical values of 
these parameter from LTE and non-LTE calculations are sum- 
marized in Lamers & Cassinelli (1999, chapter 8, "LC"). We 
have not used the absolute value of the velocity gradient in the 
line-force term, because we are interested in monotonic veloc- 
ity laws. 

The constant C represents the eigenvalue of the problem (see 
below) and is given by: 

^ , (3) 

0-EV,hM) 

where y,/, = {IkhoiT Imn)^^^ is the hydrogen thermal speed, 
is the electron number density in units of 10 '^cm^^ (Abbott 
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1982), while the meaning of all other quantities is standard 
(see, e.g., LC). 

Together with the momentum equation Q, the continuity 
equation reads: 

A-nr^p V - M (4) 

Introducing the following change of variables: 

u - -R^/r, (5) 

w = v/fl, (6) 

w' - dw/du, (7) 

where a is the isothermal sound speed, i.e., p - a^p. Replacing 
the density p from 0, the m-CAK momentum equation Q 
with the line force (|2}, we obtain: 

F{u,w,w') = il - -^\w w' + - +A (l+n^u) 

-CfDg(u){wr\wWf = (8) 



here the constants A, C and Q are: 

^ ^ GM(1 -r) ^ ^ "Ireak-up 



2n aRl 



^break-up Va' 
where is defined by: 



(9) 

(10) 
(11) 

(12) 



here v^^c is the escape velocity and vi,reak-up is the "break-up" 
velocity. The function g{u) is defined as: 



giu) = (l + Vl - m2) , 

and the constant D is: 

^ 1 {l+ZueYHe) 



(13) 



(14) 



niH (1 +AHeYHe) 

Yhc is the helium abundance relative to the hydrogen, Z/ze is the 
amount of free electrons provided by helium, Ane is the atomic 
mass number of helium and nin is the mass of the proton. 

The standard solution, from this non-linear differential 
equation (|8}, starts at the stellar surface and after crossing the 
singular point reaches infinity. At the stellar surface the differ- 
ential equation must satisfy a boundary condition, namely the 
monochromatic optical depth integral (see Kudritzki, 2002, eq. 
[48]): 



Tpho 



Jr. 



cte P (r) dr 



(15) 



A numerically equivalent boundary condition is to set the den- 
sity at the stellar surface to a specific value. 



When the singularity condition. 



-F{u, w, w') - F„' - 



(17) 



is satisfied, its location {u - Uc) corresponds to a singular (or 
critical) point, and in order to get a physical solution, the regu- 
larity condition must be imposed, namely: 

d dF dF 

—F(u, w, w')^— + —w' = F, + F„ w' = 0, (18) 
du au ow 

hereafter, all partial derivatives are written in a shorthand form, 
i.e. Fu - dF/du. 

In order to satisfy simultaneously equations (|8j, ( I17> . ( I18> 
and il5\ or il6\ . the value of the constant C is not arbitrary, i.e., 
the constant C is the eigenvalue of this non-linear problem. 



3. Coordinate Transformation and System of 
Equations 

In this section we will apply another coordinate transformation 
and introduce a new independent variable, t, without physical 
meaning. This will allow us to transform the non-linear differ- 
ential equation for the momentum into a system of coupled 
differential equations, which is numerically integrable. 
Defining 



1 2 

and 

, dy 
p = w w - — , 
du 



(19) 



(20) 



the momentum equation can be written in a general form 
as: 

F(u,y,p)^0. (21) 

Differentiating this function and using dy - p du, we obtain: 

dF ^{Fu +p F,}j du + Fpdp^O (22) 

We introduce now a new independent variable, t, defined im- 
phcitely by: 



du - Fp dr. 



(23) 



Because u and t are independent variables, they have to be re- 
lated between them. While Fp + 0, we can write t as function 
of u as: 



dy dy j du 
du dr \ dr 



= P 



(24) 



We can transform from ll22> to the following system of ordinary 
differential equations: 



p(R») -pt. 



(16) 



(25) 
(26) 
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dr 



■[Fu +pF,). 



(27) 



A solution of this system of differential equations is also a so- 
lution of the original momentum equation, since if any initial 
condition {ucyo, po) satisfies F{u^),y^), p^)) - 0, then a solution 
of (|BH57} verifies that F(u(T),y(T), p(t)) = 0. 
An advantage of this equation system ( I25H27^ over the CAK 
momentum differential equation is that all the derivatives 
are explicitely given, therefore there is no need to use root- 
finding algorithms to find the value of the velocity-gradient. 
Also, standard numerical methods (e.g., Runge-Kutta) can be 
used to integrate this system. 

4. Linearization and Eigenvalue Criteria 

All critical points of the system ( I25l27t satisfy simultaneously 
F - 0, U - and P - 0. Thus, in order to study the be- 
havior of the solution in the neighborhood of a singular point 
we linearise this system of differential equations, using the 
Groebman-Hartman theorem ("GH", see appendix. For more 
details see Amann 1990), we obtain: 

dU = (t/„ + pUy) du + Up dp, (28) 

dP = (P„ + pP,^ du + Ppdp. (29) 

We have not included dY because U and Y are dependent {Y - 
pU, eq. I26> . The GH theorem indicates that the eigenvalues 
(and the eigenvectors) of the partial derivative matrix B, defined 
by: 



B 



Uu+pUy Up 

Pu+pPy Pp 



(30) 



provide the information concerning the topology structure of 
the critical (singular) points. 

On the other hand, if we consider the eigenvalue C as a free 
paramenter, the critical points are the solution of the following 
system of equations: 

F[u,y,p,C) ^0, (31) 

U[u,y,p,c) =0, (32) 

P[u,y,p,c) =0. (33) 

In this case the number of incognits is greater than the number 
of equations, then we can only solve for the incognits in terms 
of one of them (implicit function theorem). 

If {uc,yc,Pc,C^ satisfies the previous system and further- 
more, A at the singular point, we can solve for yc, Pc and 
C and its derivatives as a function of Uc- We obtain for the gra- 
dient of C: 



dC 

duc 



-— det 

A 



U, 

Pr 



Ur. 



Fu 

p„ 



(34) 



where the determinant A is defined by: 



A = det 



(Fy Fp 
U, Up Uc 



Fc^ 



P. Pn P. 



(35) 



Considering that at the critical point: Fp - and f „ 
-pFy, equation (I34> transforms to: 



dC 1 , 
det 

dUf A 



(Fy -pFy 

Uy Up Uu 

Py Pp Pu 



Fy 



det(B) 



(36) 



The GH theorem establishes that the singular point topology is 
determined by the sign of the eigenvalues of the matrix B. A 
critical point is of X-type (Saddle) when the eigenvalues are 
both real and have opposite sign, or equivalently, 



det(B) < 0. 



(37) 



Applying this to the rotating-CAK model (see next section for 
details) it is verified that, 



> 



(38) 



then a X-type singular point corresponds to the condition: 
dC 



dUr 



< 0. 



(39) 



Therefore any X-type physically relevant singular point, is re- 
lated to the behavior of the eigenvalue C, specifically when 
condition J39> holds. 

5. The Topology of the Rotating-CAK Model 

The analysis we have shown in the last section (up to ea.l37> is 
valid for the general radiation driven stellar winds theory, i.e., 
including the finite-disk correction factor. In the remainder of 
this paper we will study the original CAK model, i.e. /d = 1. 
The topological study of the m-CAK model will be the scope 
of a future article. 

The non-linear momentum equation (|8} for the rotational- 
CAK model (including 6), in (m, y, p) coordinates reads: 

F{u,y,p) = - ^p + - + a[\ +Q}u) 

-Cg{u)y-'l^p" ^Q. (40) 

From this equation, after a straightforward calculation, the cor- 
responding system of differential equations U and P, (eg s. 1251 
and|^ respectively) are: 



1 



-(5/2 „-l+a- 



U {u,y,p) = 1 - T aCg{u)y p 



(41) 



P{u,y,p) = 



2 , p^ 
2y^ 



Cg(u)y-'-''^p'-'' 



+Cg(u)h(u) y-'l^ p" 



where 



c / 



h{u) = -6u[\-u^ + Vl -M^j , 



(42) 



(43) 
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Fy and A are: 



Table 1. B2 V Stellar Parameters 



Ly^ 2 



A = {\-a)g{u)y-'l^p" X 

((4a + ^)//2 +2Q'5Cg(M)y/?")]. 



(44) 



(45) 



It is easy to verify that both, Fy > and A > 0. 

Once the location of the critical point, m^, is known, we can 
solve yc, Pc and C from equations ( I40ll41ll42li . obtaining: 



1 a I (2 / , a 

2 (l-a)qc\Uc ^ 'j 

Pc = \qc + ( — + A (l + Q^,.) I , 

C = — +a{\ +D}u,] \ , 

where (7^. is the positive solution of the quadratic equation 

+ 7-^ ( — + A (1 + Q-M,.) I ^ = y (mc) , 
with 



(47) 



(49) 



r(Mc) = -2AQ^ + — 



2/z(mc) 1 2 / ^2 \\ 
^--.A(l+Q«,)j. 



(50) 



Equations j46ll47l andl48> are generalizations of the equations 
[49], [50] and[51] from Marlborough & Zamir (1984), includ- 
ing now the effect of the line-force term (ueIWY . 
Since the eigenvalue C must be positive, we have: 



'2 

— +a(1 -hQ^Mc) > 0. 



(51) 



This last inequality imposes a restriction in the position (fur- 
thest from the stellar surface in the radial coordinate r) of the 
singular point, namely: 



(52) 



a(i + Vi -8oVa) 
6. Topology of the rotating CAK wind 

In this section we show the results of our topological analysis. 
Following Bjorkman (1995), we choose a typical B2 V star with 
stellar parameters summarized in tabled We have adopted the 
line-force multiplier parameters k and a from Abbott (1982) 
and the S parameter is from LC, these values are summarized 
in table 12 



RIRq 


M/Mq 


L/Lq T.fflK r A 


4.5 


9.0 


3553. 21000. 9.27 10-3 1413. 



Table 2. Line-force Parameters 



k a 6 
0.212 0.56 0.02 



(46) 6.1. The frozen-in ionization (S -0) 



The factor (riElWf in the line-force takes into account the 
changes in the ionization of the wind. As a first step to under- 
stand the topology of the rotating-CAK model we set 5 = 0. 



(48) 6.1 .1 . The critical point interval 

This case is simpler because we can obtain analytically from 
equations J46I47I481 the variable y as function of u and p, i.e.: 



2 (2/m + A{\+a^u)+p-Cp"')' 



(53) 



Bjorkman (1995) obtained the same result (see his equation 
[13]) for the non-rotational case, Q = 0. Moreover, qc can be 
expressed by: 



qc 



j{uc) = 4m-^ - 2A Q- 



(54) 



As we mentioned in previous sections, qc must be positive. This 
gives an extra restriction for the location of the singular point 
(nearest to the stellar surface in the radial coordinate r). The 
value of M„„„ is defined when j{umin) - and is given by: 



max 



-1, 



/27a 



(55) 



Therefore, Uc is restricted to the interval (m„„„, «,„«.«:)■ A similar 
result which restricts the location of the critical point to an 
interval has been found by Marlborough & Zamir (1984) see 
their eq. [65]). 

Figure ^ shows the behavior of y (eq. I54> versus rjR^ - 1 
for different values of the rotational parameter Q.. For the 
non-rotational case Q = 0, the function y is positive for 
almost the whole integration domain, therefore the singular 
point can be placed anywhere. Thus, is the lower boundary 
condition which fixes the position of the singular point. For 
the rotational case, the second term in the RHS of equation 
(I54> is the dominant term for almost any value of Q. Therefore 
the larger is Q the larger is the value of r^m (= -^/umin) as the 
diff'erent curves in figure^ show. 

On the other hand, the value of the term 8Q^/A in eq. i52\ 
is almost negligible and consequently the value of rmax 
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Table 3. Analytical approximation and numerical results for 
the rotational-CAK model with 6-0. Note: the mass loss rate 
is given in units of 10"^ MQ/year and the terminal velocity is 
in km/ sec 



Q 


r,„m/R, 


rmaxIR. 


rcritlR. 




C 


M 







1 


706.50 


1.559 


48.248 


48.223 


1.814 


1008.2 


0.2 


5.316 


706.46 


5.370 


47.972 


47.961 


1.831 


920.1 


0.4 


10.632 


706.34 


10.669 


47.648 


47.634 


1.854 


800.3 


0.6 


15.948 


706.14 


15.992 


47.322 


47.303 


1.877 


683.3 


0.8 


21.264 


705.86 


21.334 


46.992 


46.964 


1.902 


574.5 



jcr 




f R. -J 



Fig. 1. a) The function y versus r//?, - 1 ('ea.l54> for different 
values of the rotational parameter Q. O = 0, continuous-line; 
Q - 0.25, dashed-line; Q = 0.8, dotted-line. The horizontal 
line is for y = 0. 

b) The eigenvalue from eq. ( I48> as a function of r/R^ - 1 for 
the same values of Q as a). Note that dC/du is always negative 
in the interval (r„„„, r„,a.f). 



(= -\ I Umax) is almost constant at r = AR^/2. Table |3l shows 
the interval r„,ax) for different values of the parameter 
Q.. It is clear from this table that from a very low value of 
f2 - V2/A (= 0.038 for our test star) the location of r„„„ > /?» 
and the critical point, rent, is strongly shifted downstream in 
the wind. Once r„,Y (or Uc) is fixed, the value of yircru) is 
inserted in eq. \5A\ and qc and C are obtained. FiguredJi shows 
C from eq. J48I I against - 1 for the same values of Q. of 
figure^ (same type of lines too). The value of the derivative, 
dCldu is always negative indicating that the critical point is an 
X-type. 



6.1.2. Linearization of tine B-IVIatrix: Eigenvalues and 
Eigenvectors 

The matrix of linearization B (eq.13011 is given by: 



B\ 



q2 a(\-a)ulri ^ 



(56) 



where 6 and rf are given by: 

6»= ^l-Q?Au^J2, (57) 
and 

77 = -(2+Am,(1+q2m,))/(1-Q') ■ (58) 
The eigenvalues of the Z?-matrix are: 

yu± = + J9&^+4-ail-a)(l+20^)j]] . (59) 

It is easy to check that the product of the eigenvalues fi+ x yu_ 
is negative, i.e., the topology of the singular point is saddle 
or X-type as the GH theorem establishes. We can write the 
associated eigenvectors as (1, v±), where 



a (I - a) u^i] 



(60) 



where v+ (v_) corresponds to the unstable (stable) manifold. 
The maximum value of C, that accounts for the minimum mass 
loss rate, occurs when u^ - Wmin, eq- (I48> becomes: 



A(i-n -)- 



+ j^{A(l-Q?)-2y" , ifD. < V27A 



(61) 



6.1 .3. Pliase Diagram 

In order to obtain a numerical solution for the wind, we need 
an approximation for the location of the critical point. We use 
the value of u^ - u„i„ as a first guess of the critical point and 
use this value of u^ to calculate C„„„ as our first guest for the 
eigenvalue C. Table|3shows the values of r^m, ?"«■,( and C,„,„, C 
confirming that this is a very good approximation. Now using 
standard numerical algorithms we integrate our system: U - 
0, y = and P = (eqs. |25]|26land|27] respectively) from 
the singular point up and downstream to obtain the numerical 
solution. 

As Bjorkman (1995) pointed out, it is insightful! to study 
the solution topology in a p versus (r//?, - 1) diagram. Figure 
IJlshow this phase diagram. If we start to integrate at the singu- 
lar point, we cannot leave this point because all the equations, 
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r R, -J. 



Fig. 2. The topology of the freeze in ionization rotating-CAK 
model (6 = 0), p versus {r/R, - 1) for Q = 0.25. The unique 
curve that starts at the stellar surface and reaches infinity is the 
CAK original solution (continuous-line). 

U - 0, Y - and P = are simultaneously satisfied at this 
point. Therefore we have to move slightly up and downstream 
along the direction of the unstable manifold. After this, we can 
integrate obtaining the different solutions showed in figure |2] 
From this figure, it is clear that the only solution that reaches 
the stellar surface (r — > +oo, in the direction of -(1, v+)) and 
also reaches infinite (t — > +oo, in the direction of (l,v+)) is 
the original CAK solution (continuous-line). The results from 
table |3]for the non-rotational case are the same one obtained 
by Bjorkman(1995). 




in- J irff j£J" JO* 



Fig. 3. The velocity profile v(r) (in km/ sec) as function of 
r/Rt - 1 for the non rotating case (continuous-line) and for the 
rotational cases: Q = 0.25 dotted-line; Q = 0.8 dashed-line. 
The location of the singular points is indicated by a larger-dot 

Figure |3l shows the velocity profile, v (in km/ sec) versus 
r/Rt - 1 for our B2 V test star We have chosen this value of 
= 0.25 from the study of Abt et al. (2002), that concluded 
that B-stars rotate at a 25% of their break-up speed. The 
values Q. = 0.8 accounts for a fast rotator, e.g., a typical 
Be-Star (Chauville et al. 2001). We see from this figure that 
neglects the rotational speed always overpredicts the value of 
the terminal velocity. 



We conclude that the rotational speed shifts the location of 
the critical point downstream and reduces the terminal veloc- 
ity, but has almost no influence on the value of the eigenvalue 
(mass loss rate). Furthermore, we can see from our approxi- 
mate and numerical results summarized in tablel^that the CAK 
wind do not have the smallest possible eigenvalue, C„„„, or the 
maximum mass-los rate as Feldmeier et al. (2002) and Owocki 
& ud-Doula (2004) concluded for a non-rotating CAK model 
with zero sound speed. Contrary to expectation, the rotating- 
CAK wind critical solution corresponds to an almost minimum 
mass-loss rate {maximum eigenvalue). 

6.2. The rotating CAK model (6 i= 0) 

Abbott (1982) studied the indirect influence of the density 
in the line-force through the dependence of the ionization 
balance in the electron density. He found that this dependency 
modifies the force-multiplier and therefore the line-force by 
a factor (ns/W)*, where is the electron density (in units of 
10"" gr /cm^) and W is the dilution factor. 
Although this is a weak influence, because 6 ranges between 
0.0 and 0.2, it is important to study how its inclusion in the 
momentum equation modifies the topology of the rotating 
CAK model. 

As we pointed out in section |5] the existence of q,. > 
implies y (Uc) > 0. The behavior of y (eq.|50j versus r/R, - 1 
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Fig. 4. Same as figure^a) and b) but for 6 - 0.02. See text for 
discussion. 
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is shown in figure|3t for the same values of Q. of figure[2 We 
can interpret from this figure, that y in (R,, rmax) is a decreasing 
function in the neighborhood of r ^ and an increasing 
function in the neighborhood of r ^ r^ax- Furthermore, 
the function 7 posses only one minimum in the integration 
domain located at r = rmin- The value of y{r„in) decreases as Q 
increases. For small values of Q, y(r„,j„) > and the location 
of the singular point rent can be anywhere in the integration 
domain {Rt,r,„ax)- For larger values of y(rmm) < and 
y(r) = has two roots, at r = r± where r_ < r+ (or u- < u+). 
Thus, rcrii can be located now in two different intervals, i.e., 
Kcrit G {Rt,rJ) or Kcrit 6 ('"+, '"ma.i)- For the particular case 
where yir„i„) = 0, the value of the rotational speed parameter 
is Q = Q.bif, where the subscript 'bif accounts for the 
bifurcation in the wind solution topology. 

Figure]^ shows C (eg. I48l l against r/Rt - 1 for different 
values of Q.. When y(r) is positive, C exhibits now a different 
behavior compared with the 6-0 case. As long as Q < Q^,/, 
dC/du is positive in the interval {Rt,r,„i„), i.e., any singular 
point in this interval is an attractor. Furthermore, dC/du 
reaches its maximum, when y is minimum, i.e., in the neigh- 
borhood of r = r„„„, from this point up to r,„ax, dC/du < and 
any singular point in (r„j„!, rmax) is X-type. When Q > Q.bif, the 
minimum of y{r) is negative and the intervals are reduced to: 
(Rt,r-) for the attractor type singular point and {r+,rmax) for 
the X-type singular point. 



The behavior of jc, pc and C (eas. 1461 1471 and l48t in the 
neighborhood of r — > r±, in the inverse radial coordinate m, is 
as follows: 



Table 4. Numerical results for the rotation-CAK model with 
6 - 0.02. Note: the mass loss rate is given in units of 
10"^ MQlyear and the terminal velocity is in km/ sec 



Q 




C 


M 







4.712 


50.949 


1.644 


929.6 


0.2 


7.196 


50.726 


1.657 


879.3 


0.4 


11.926 


50.291 


1.683 


779.0 


0.6 


17.145 


49.790 


1.713 


669.9 


0.8 


22.576 


49.247 


1.747 


564.6 



We conclude from that the factor (riE/Wy strongly shifts 
outwards the location of the critical point and produces a bifur- 
cation in the solution topology. 




yc 

Pc 

c 



\U — U±\ ' — > OO . 

a I 2 



\ - a \ u± 



H-A(l -hQV) 

+ 00 



(62) 
(63) 
(64) 



This divergent behavior of C in the neighborhood of r-t, is 
shown in the curves for Q 9^ in figure|4] The almost constant 
value of C explains why a slightly change in the eigenvalue 
cause an enormous change in the location in the singular point. 

Tables|3]and|5lsummarise the numerical calculation for our 
test star with 6 - 0.02 and 6-0.1 respectively. The data of 
the M column (see also the C column) show that the effect of 
the rotation on M (C) is almost negligible. Figure |5] shows the 
velocity profile for three different values of £1 (0.0; 0.25; 0.8), 
panel a) for 6 - 0.02 and panel b) for 6 - 0.1. A large dot 
shows the respective positions of the critical points. The effect 
of shifting the position of the critical point is stronger for 
low rotational speeds and decreases when Q increases as 
a comparison between figures |3 and |5] clearly shows. The 
terminal velocity, is a decreasing function of the rotational 
speed and has almost the same behavior as in the <5 = case. 
But for high rotational speeds, the influence of 6 in Vco is 
negligible/small as a comparison between tableland table l4l5l 
shows. 




Fig. 5. The velocity profile v(r) (in km/ sec) as function of 
r/R, - 1 for the non rotating case (continuous-line) and for the 
rotational cases: Q. = 0.25 dotted-line; Q = 0.8 dashed-Une. 



6.2.1. Bifurcation rotational speed 

In order to have an analytical approximation for the value of 
Qfc,y, we can approximate for the function h(u), eq. (I43> : 



h(u) 



-u , 



(65) 
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Table 5. Same as Table Q but for 6^0.1. 



Q 


rcritlR. 


C 


M 







11.708 


63.099 


1.122 


801.1 


0.2 


13.168 


62.734 


1.133 


775.1 


0.4 


16.873 


61.876 


1.162 


708.5 


0.6 


21.853 


60.619 


1.205 


621.1 


0.8 


27.683 


59.107 


1.261 


528.5 



then the minimum of y{u) is achieved at: 



1 



Ubif 



6A 



1/3 



(66) 



and the minimum value of y{uhif) is 

j(uhif) ^ -2 (1 - a) q2 A (l + 2 (1 - a)-^'^ d^'^A-^'^) 

+3 {\-ayi^6^I^A^I^ -25 (67) 
0, we can obtain the bifurcation value of Q. 



From yiubif) 
given by 



^bif - 



3 (1 - a) '^^(^2/3^2/3 _ 2 

^ 2 (1 - a) A (l + 2 (1 - a)-^l^ 6^I^A-^I^) 



(68) 



Figure|6lshow the curves r„„„ and r-t as function of D.. The 
intersection point for all the curves (in continuous-line) is at 
Q - Q.bif- Critical points can not be located between the curves 
r+ and r_ (filled region). In addition, we show curves for 
the location of the critical point from numerical calculations 
(dashed-Unes), with the lower boundary condition, Tp/,0, = 2/3. 




n.z n.* a. 




Fig. 6. The curves r„„>, and r± are shown with continuous lines 
as function of Q. These curves represent the boundaries for the 
type of the singular point topology. The numerical results for 
the location of the singular point is also shown in dashed-Une. 
Figure a) is for S - 0.02 and b) for 6 = 0.1. See text for details. 



We can clearly see from this figure, that the position of the 
critical point is shifted outwards from the stellar surface and 
the greater is S the further is the position of this critical point. It 
can be inferred from this figure, that the location of the singular 
point remains almost constant as long as Q < Quf, but from 
values of O > Q;„y the position of r„.„ grows almost linear 
with Q. 

This behavior of the solution topology can be applied for 
the winds of Be-Stars. At polar latitudes, i.e., slow rotational 
speed, the wind behaves as the standard CAK wind, but as the 
latitude approaches to the equator, the rotational speed is larger 
than Q.bif and the wind is slower and denser. This transition 
from polar to equatorial latitudes seems to have a similar be- 
havior described by Cure (2004) for the more general rotating 
m-CAK wind. The study of the influence of this bifurcation in 
the winds of Be stars will be the scope of a forthcoming article. 



(Feldmeier et al. 2002). As we mentioned in previous section, 
we move slightly from the singular point in the unstable mani- 
fold in the directions ±(1, v+) and then integrate outwards and 
inwards, obtaining the solution topology of figure Figure 
|8] shows both shallow and steep solutions, but here we have 
started from the singular point but with different values of the 
eigenvalue C. The almost horizontal dotted curves correspond 
to pc from eq. (147 > and is not a continuous curve because of 
the bifurcation in the solution topology. The left-most shallow 
curve that start at a X-type singular point do not reach stel- 
lar surface, but is trapped by the attractor singular point. This 
result implies that there is a shorter window of locations of sin- 
gular points, that fixes the eigenvalues C and can reach infinity 
passing through the X-type critical point. This feature is not 
present when 6 -Q. 



6.2.2. Phase Diagram with 5i^0 . 

Figure shows the phase diagram p versus r/R, - I for D. = 
0.25 and 6 = 0.02 for our test star. The solution topology seems 
to be similar to the i5 = case. Here the shallow and steep 
curves are in continuous line and dashed line respectively. The 
shallow solution is the CAK solution while the steep solution 
coiTespond to accretion flows or for radiation driven disk winds 



7. Conclusions 

In this article we have examined the topology of the rotating- 
CAK wind non-linear differential equation. After the introduc- 
tion of an additional physically meaningless independent vari- 
able (t) we transform the momentum equation to a set of equa- 
tions where all the derivatives are explicitly given. This for- 
malism permitted us to linearise the equations in the neighbor- 
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Fig. 7. The topology of the rotating-CAK model, p versus 
r/Rt - 1 for Q = 0.25 and 6 - 0.02. The unique curve that 
starting at the stellar surface and reaches infinity is the CAK 
original solution (continuous-line). 




Fig. 8. The topology of the rotating-CAK model, p versus 
(r/Rt - 1) for Q = 0.25. The different curves correspond to 
different values of the eigenvalues. The left-most curve starts 
at the 'normal' critical point and is trapped by the attractor See 
text for details. 

hood of the critical points. This linearization let us to define a 
condition for the derivative of the eigenvalue that defines the 
topology of the critical point, i.e, X-type or attractor. 

We have applied our results to the case of a point star 
(CAK) for a frozen in ionization rotating wind, recovering and 
generalizing the results of previous studies (Bjorkman 1995 
and Marlborough & Zamir 1984). The most significant result 
(with 5 = 0) is that the wind does not assume the maximum 
mass-loss rate but almost the minimum. 

For the more general case, where changes in the wind ion- 
ization are taken into account, our analysis shows the existence 
of a bifurcation in the solution topology, where two critical 
points exist. The first critical point (closer to the star's surface) 
is an attractor while the second is the standard CAK critical 
point. Besides the known fact that the rotational speed shifts the 
location of the critical point outwards in the wind, the inclusion 
of the term (n^/W)'' produces the same effect, reinforcing this 
displacement. 

The bifurcation topology seems to explain the results from 
Cure (2004) that there exist two regions in the wind of a fast 



rotating hot star, one where the wind is the one from the stan- 
dard solution (fast wind) and the other with a new solution that 
is slower and denser This result shows us the necessity to per- 
form a topological analysis of the rotating m-CAK wind. This 
study is currently underway. 

Appendix A: Linearization in tlie neighborhood of 
a singular point 

Considering the system of differential equations given by: 



y 



dx/dt 
■■ dy/dt 



X(x,y), 
--Y(x,y), 



(A.l) 



where x and y are the dependent variables and t is the indepen- 
dent variable. The point (xcje) is singular (critical) if and only 
if verifies that: 

Xix„y,) =0 
Y(x„y,) =0 

In order to understand the topological behavior of the solutions 
close to a singular point, we expand this system using Taylor 
series at (xc,yc), obtaining: 

■* = + ^x\ix„y,) (x - Xc) + Xy\(j,^,y^) (y -yc) + o (x,y) 

y = J'kx^y,) + Yx\(x,.y,) (X - Xc) + Yy\^x„y,) (j - Jc) + O (x, y) 

Neglecting superior order terms and using that (xc,yc) is a sin- 
gular point, we have, in matrix form: 



(A.2) 



B = 



B 

y I \y-yc 

where B is the Jacobian matrix at (xcje), i.e., 
Xx X, \ 

For a 3-dimensional system case: 

X - X(x,y,z) , 
y = Y(x,y,z) , 
z -Z{x, y, z) , 

with a constraint (p {x, y, z) satisfying: 

(l)xX + 4>yY + 4>,Z = , 

the surface defined by: 

S ^{(x,y,z):<f>(x,y,z)^0}, 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



is invariant for the differential equation system evolution. 

From the equation (IA.6> . we can solve for z (using the im- 
plicit function theorem) obtaining z - z(x, y) and reduce the 
system (IA.4> to: 



X = X{x,y,z (x, y)) , 
y = Y(x,y,z(x,y)) . 



(A.7) 
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If (xc, yc, Zc) is a singular point of (IA.7> . the Jacobian matrix 
B is given by: 



B = 



Y, + z,Y, Yy+ZyY, 
Since Zx - -4>x/4>z and Zy 



(Av,y,) 



-(py/(pz, we have 



1 

5 = — 

4>z 



ixY, <p,Yy 



byY, 



(A.8) 



(A.9) 



{xc.yc) 



For <p - F, X = U and Y - P,we obtain equation (I30> . 
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